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Ringdal's  maximum  likelihood  method  of  removing  magnitude 
bias  was  tested  by  removing  the  apparent  bias  of  surface  wave  magnitude 
estimates.  Bias  removal  was  convincingly  demonstrated  by  comparing  maxi- 
mum likelihood  estimates  of  obtained  by  a single  sensor  to  those  obtained 
by  an  array  at  the  Alaskan  Long  Period  Array  (ALPA)  site.  Since  the  beam- 
formed  array  has  a lower  detection  threshold  than  the  single-sensor  reference 
site,  it  can  serve  as  the  standard  by  which  to  judge  whether  positive  magnitude 
bias  has  been  removed  from  the  reference  site  surface  wave  magnitude  esti- 
mates. 

The  results  of  this  test  indicate  that  below  the  90  percent  detec- 
tion threshold  of  the  reference  site,  positive  magnitude  bias  becomes  a sig- 
nificant factor  in  surface  wave  magnitude  estimates.  Comparison  with  beam- 
formed  array  data  indicates  that  application  of  Ringdal's  method  accurately 
removes  this  bias  down  to  approximately  the  reference  site  25  percent  detec- 
tion threshold.  Below  this  point,  the  method  over- corrects  the  magnitude 
estimates,  resulting  in  abnormally  low  values.  This  test  also  indicates  that 
correcting  surface  wave  magnitude  estimates  for  bias  by  Ringdal's  method 
will  improve  the  value  of  the  Mg-  discriminant  for  network  data  by  making 
it  applicable  to  events  with  lower  bodywave  magnitudes.  The  results  obtained 
in  deriving  an  unbiased  Mg-  m^  relationship  were  consistent  with  various 
other  published  results  and  explain  their  diversity  by  showing  that  the  slope 
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SECTION  I 
INTRODUCTION 


A.  THE  MAGNITUDE  BIAS  PROBLEM 

The  conventional  method  of  estimating  seismic  event  magni- 
tudes is  to  average  the  magnitude  measurements  made  at  those  stations  of  a 
network  which  detect  the  event.  The  purpose  of  this  averaging  process  is  to 
'smooth  over1  any  abnormally  low  or  high  deviations  of  the  magnitude  esti- 
mates caused  by  source  effects  such  as  event  radiation  pattern,  regional 
magnitude  effects  due  to  upper  mantle  structure  under  the  station  (e.  g.  , Der, 
1977),  local  receiver  site  geology  and  crustal  structure  effects  (e.  g.  , Ring- 
dal  et  al.  , 1973;  and  Chang  and  von  Seggern,  1978),  and  abnormal  path  at- 
tenuation effects  (e.  g.  , Der,  1977;  and  Sun,  1977).  However,  when  estimat- 
ing the  magnitude  of  events  near  the  network  detection  threshold,  the  averag- 
ing process  itself  is  a biased  procedure  due  to  seismic  noise,  which  limits 
signal  detection  to  those  few  stations  of  the  network  which  have  the  largest 
positive  amplitude  deviations. 

It  is  important  to  note  that  the  magnitude  bias  to  be  discussed 
in  this  report  is  not  due  to  the  interference  of  noise  with  weak  signals.  It  is 
assumed  that  detected  signals  are  measured  at  high  enough  signal-to-noise 
ratio  to  render  negligible  the  influence  of  noise  interference  on  magnitude 
measurement.  The  magnitude  bias  discussed  here  is  grounded  in  two  other 
quite  different  physical  phenomena  - the  effect  on  signal  detection  of  source- 
to- station  statistical  variation  in  signal  amplitude  for  seismic  events  due  to 
source,  path,  and  local  receiver  effects,  and  the  temporal  statistical  variation 
of  the  seismic  noise  levels  at  the  stations.  These  factors  influence  the  signal- 
to-noise  ratio  variation  between  observations  which  in  turn  determines  that  only 
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the  larger  signals  from  an  event  are  used  to  estimate  the  event  magnitude. 
The  net  result  is  a positive  bias  in  the  event  magnitude,  which  increases  as 
the  true  event  magnitude  decreases. 


- • 


A method  for  removing  this  positive  magnitude  bias  was  de- 
veloped by  Ringdal  (1975).  In  brief,  his  method  is  to  measure  signal  magni- 
tudes at  the  detecting  stations  and  to  determine  upper  bounds  on  the  signal 
magnitude  at  the  non-detecting  stations  by  measuring  the  magnitude  of  noise 
at  the  station  as  though  it  were  a signal.  A maximum  likelihood  estimate  of 
the  true  network  magnitude  and  standard  deviation  is  then  computed  from  the 
combined  measurements  of  noise  at  non-detecting  stations  and  measured  sig- 
nal magnitudes  at  detecting  stations. 

Any  test  of  Ringdal1  s method  must  of  necessity  establish  that 
bias  in  the  measurement  of  the  magnitude  of  small  signals  exists  as  a signif- 
icantly observable  effect  and  can  be  effectively  removed  down  to  some  mini- 
mum level  of  detection  probability.  In  light  of  this,  the  results  of  a paper  by 
Evernden  and  Kohler  (1976)  should  be  mentioned.  The  authors  of  this  paper 
investigated  the  question  of  bias  in  estimates  of  bodywave  magnitude  at  small 
magnitudes.  Their  results  indicate  that 

• Network  bodywave  magnitudes  do  have  the  potential  of  being 
biased. 

• The  probability  of  any  significant  bias  occurring  is  so  small 
that  it  can  have  only  negligible  effects  on  use  of  an  Mg-  m^ 
discriminant. 

It  therefore  appears  more  profitable  to  use  surface  wave  magni- 
tudes rather  than  bodywave  magnitudes  in  this  study  of  positive  magnitude  bias. 


B. 


THE  TASK 


The  primary  goal  of  this  report  is  to  test  Ringdal's  method  of 
correcting  surface  wave  magnitude  bias.  In  the  course  of  carrying  out  this 
test,  the  following  points  will  be  covered: 

• Design  an  experiment  to  carry  out  this  test 

• Determine  the  differences  between  conventional  network  magni- 
tude estimates  and  the  estimates  derived  from  Ringdal's  meth- 


• Describe  qualitatively  the  conditions  under  which  Ringdal's 
method  can  be  expected  to  function  properly 

• Determine  the  effect  of  Ringdal's  method  on  the  M - m,  dis- 

s b 

criminant. 

Section  II  of  this  report  discusses  the  methodology  of  Ringdal's 
method  of  magnitude  estimation,  the  experiment  designed  to  test  this  method, 
and  the  data  base  used  to  carry  out  the  experiment.  Section  HI  presents  an 
analysis  of  the  results  of  performing  the  experiment  in  terms  of  the  Mg- 
relationship  and  the  m^  discriminant.  Section  IV  presents  the  conclusions 

of  this  report.  Section  V lists  the  references  cited  in  this  report. 


SECTION  II 
METHODOLOGY 
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A.  RINGDAL'S  NETWORK  MAGNITUDE  ESTIMATION  METHOD 

Consider  a large  seismic  event  detected  by  each  of  the  stations 
of  a network.  At  each  station,  the  magnitude  of  the  event  is  measured.  From 
Freedman  (1967),  one  would  expect  these  station  magnitudes  to  be  normally 
distributed.  Therefore,  with  an  unbiased  selection  of  a large  enough  sample 
of  magnitude  measurements,  the  average  of  these  station  magnitudes,  called 
the  network  magnitude  estimate,  would  accurately  represent  the  true  magni- 
tude of  the  event.  However,  as  the  true  magnitude  decreases  toward  the  net- 
work detection  threshold,  those  signals  representing  the  lower  side  of  the 
normal  distribution  are  not  detected  due  to  relatively  higher  seismic  noise. 
Averaging  the  magnitudes  of  the  remaining  portion  of  the  distribution  results 
in  a positive  bias  in  the  network  magnitude  estimate  relative  to  the  true  magni- 
tude of  the  event. 

Ringdal's  method  for  removing  positive  bias  in  the  measure- 
ment of  the  magnitude  of  a seismic  event  is  to  measure  the  magnitude  of  the 
propagated  signals  at  all  detecting  stations  and  estimate  an  upper  bound  on 
the  magnitude  of  signals  at  all  non-detecting  stations.  These  upper  bounds 
are  obtained  by  measuring  the  largest  noise  peak  in  the  signal  window  and 
computing  the  corresponding  'noise  magnitude'.  The  statistical  method  of 
maximum  likelihood  estimation  is  then  applied  to  these  signal  and  noise  magni- 
tudes to  estimate  the  unbiased  network  magnitude. 

This  method  takes  as  its  basic  assumption  that,  for  a given 
event,  world-wide  magnitudes  follow  a Gaussian  distribution  with  mean  fi 
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and  variance  O . The  unknown  mean  of  this  distribution  is  called  the  true 
network  magnitude  of  the  event.  It  is  also  assumed  that,  at  a given  station, 
an  event  is  detected  if  the  station  magnitude  m exceeds  a certain  threshold 
value  a . Thus,  the  probability  of  detecting  a given  event  may  be  written  as 


where 


and 


P(detect/ H,o)  - P(m>a/M»CT)  = 4^—^ 
= f <t>( t)  dt 


4 


,,  , 1 . -x  /2 

<p(x)  = — 7=  * e 


y/z IT 


(n-i) 


(n-2) 


(n-3) 


Now  assume  that,  for  a given  event,  records  from  an  n- station 
network  are  examined.  Further,  assume  that  the  detection  thresholds  a., 
i = 1 , 2,  . . . , n are  known,  and  that  a magnitude  m^  is  computed  at  each  of 
the  detecting  stations.  Finally,  assume  that  all  station  observations  may  be 
considered  to  be  independent. 

Given  this  set  of  assumptions,  one  computes  as  a function  of 
the  unknown  event  parameters  (H,0)  the  likelihood  function  for  a given  set  of 
signal  and  noise  magnitudes  using  as  starting  points  a preset  value  of  O and 
the  mean  signal  magnitude  for  fi  and  then  maximizes  this  likelihood  function 
with  respect  to  the  unknown  parameters.  The  likelihood  function  is: 


L(mlt . . . = 


n [-f  *+(-v-)]  * n *(-ir)  <“-4> 

all  L \ /J  all  non.  \ / 


detections 


all  non- 
detections 


where  ^ and  4 are  as  defined  in  equations  (II-2)  and  (II-3).  Note  that  the 
first  group  of  products  of  equation  (II-4)  represents  the  likelihood  function 
when  all  non-detections  are  ignored,  and  is  maximized  by  the  conventional 
estimate  of  H , the  average  of  the  magnitudes  from  the  detecting  stations. 
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The  procedure  is  then  to  maximize  equation  (II- 4 ) in  terms  of 
the  unknowns  /i  and  a . This  is  carried  out  by  an  iterative  procedure  in  the 
computer  program.  Fortunately,  it  appears  the  parameter  a can  be  re- 
stricted to  within  a fairly  narrow  range.  For  example,  Vieth  and  Clawson 
(1972)  found  a value  of  a = 0.  4 to  be  representative  of  the  WWSSN  network. 
Bungum  and  Husebye  (1974)  showed  that  a value  of  a = 0.  3 appeared  to  be 
independent  of  event  magnitude. 

A formal  development  of  this  method  may  be  found  in  Ringdal' s 
report  (1975).  Ringdal  demonstrated  that  his  method  effectively  removed  bias 
in  estimates  of  network  m^  by  simulating  normal  magnitude  errors  and  noise. 
Ringdal' s method  presumes  that  the  magnitude  deviation  of  signals  are  inde- 
pendent of  the  noise  at  the  seismic  stations  where  the  signals  are  observed. 
There  is  some  observational  evidence  that  there  is  a strong  dependence  on 
regional  magnitude  deviations  of  signals  and  noise  at  the  seismic  stations  as 
shown  by  Dietz  and  Sax  (1977)  which  leads  to  nearly  negligible  m^  bias;  but 
that  random  variations  due  to  local  geology  and  source  factors  are  probably 
much  larger  as  shown  by  Chang  and  von  Seggern  (1978).  It  is  assumed  here 
that  no  such  complication  exists  in  using  Ringdal 's  method  to  remove  bias  in 
estimates  of  long-period  M . 

B.  THE  EXPERIMENTAL  PLAN 

In  designing  an  experiment  to  test  Ringdal's  network  magnitude 
bias  correction  method,  one  immediately  encounters  the  problem  of  deter- 
mining a standard  by  which  to  judge  whether  the  bias  has  been  removed.  As 
originally  formulated  by  Ringdal,  this  requires  a perfect  network  (all  stations 
detecting  all  events  down  to  the  network  detection  threshold)  to  provide  the 
standard  unbiased  network  magnitudes.  Since  no  such  network  exists,  the 
experiment  was  rearranged  by  interchanging  the  roles  of  stations  and  events. 
All  observations  of  long-period  signals  were  made  at  the  Alaskan  Long  Period 
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Array  (ALPA).  The  assumption  that  averaging  magnitudes  of  long-period  sig- 
nals is  the  major  source  of  bias  effects  in  the  relationship  rather  than 

bodywave  magnitudes  was  prompted  by  the  observations  of  Section  I,  which  in- 
dicate that  magnitude  bias  may  be  more  evident  in  surface  wave  magnitudes 
than  in  bodywave  magnitudes.  The  approach  will  therefore  be  to  use  Ringdal's 
method  on  ALPA  data  to  estimate  the  average  surface  wave  magnitude  of  many 
events  having  the  same  bodywave  magnitude.  Since  the  beamforming  process 
reduces  the  50  percent  detection  threshold  of  ALPA  by  approximately  one-half 
m^  unit,  for  that  range  of  m^  the  beamformed  array  data  will  detect  essentially 
all  events  while  the  single- sensor  reference  site  of  the  ALPA  array  will  detect 
less  than  half  of  the  events.  Thus,  for  this  range  of  m^  values,  the  array  can 
serve  as  the  standard  against  which  to  judge  whether  Ringdal's  method  has  re- 
moved the  magnitude  bias.  (Since  the  beamforming  process  attenuates  the  sig- 
nal slightly,  one  can  expect  the  reference  site  unbiased  magnitudes  to  be  slight- 
ly greater  than  the  beamformed  array  magnitudes.  ) 

This  approach  of  course  requires  that  Mg  values  recorded  at  a 
given  station  for  events  of  the  same  m^  be  normally  distributed  as  are  the  m^ 
values  of  a given  event  recorded  by  a network.  To  examine  this,  event  minus 
mean  signal  Mg  differences  were  computed  for  all  events  above  the  90  percent 
detection  threshold.  These  values  were  normalized  by  the  transformation 


yi  = 


x.  - x 
1 


(II- 5) 


where 


x^  = event  for  a given  m^ 
x = mean  signal  at  a given  m^ 

° = standard  deviation  associated  with  x . 

The  y.  values  were  grouped  into  0.  2 unit  bins,  producing  the  histogram  in 
Figure  11- 1 where  the  data  were  normalized  so  that  the  highest  bar  had  a 
value  of  one.  A normal  distribution  (zero  mean,  unit  standard  deviation)  is 
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also  shown  in  this  figure.  This  result  indicates  only  small  deviations  between 
the  observed  histogram  and  the  presumed  normal  distribution. 

C.  THE  DATA  BASE 

The  source  of  the  data  used  to  test  Ringdal's  magnitude  bias 
correction  method  is  a suite  of  events  processed  and  analyzed  during  the 
evaluation  of  the  Alaskan  Long  Period  Array  (Strauss,  1973).  A total  of  233 
events  with  epicenters  in  the  Kurile  Islands -Kamchatka  area  were  selected 
from  the  available  ALPA  data  base.  Events  from  these  regions  which  were 
not  selected  were  rejected  due  to  the  presence  of  interfering  waveforms. 

Plots  of  these  events  were  examined  to  determine  which  events 
had  detectable  surface  waves  at  ALPA.  Surface  wave  magnitudes  were  com- 
puted for  each  event,  whether  detected  or  not,  by  a program  which  searches 
the  signal  gate  for  zero  crossings  and  the  interval  between  each  pair  of  zero 
crossings  for  the  maximum  absolute  amplitude.  Using  these  amplitudes  and 

the  period  computed  from  the  zero  crossings,  a series  of  M values  was  com- 

8 

puted.  This  series  of  values  was  then  searched  for  the  maximum  20  second 

period  M values.  This  method  of  computing  M is  described  in  detail  in  an 

8 8 

earlier  report  (Strauss  and  Weltman,  1977).  The  only  change  made  is  that, 
following  Nuttli  and  Kim  (1975),  the  Mg  values  used  in  this  report  were  com- 
puted from 

Ms  = log10  Y + l'07  l°810  A+  l-  16  10°^  A <25°  (II-5) 
Ms  = logio  4"  + 1,66  logio  A+  °* 30  A>25°  (n-6) 

where 

A a peak-to-peak  amplitude  in  millimicrons 
T = period  in  seconds  measured  between  17  and  23  seconds 
A a epicentral  distance  in  degrees. 
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In  order  to  'bin'  the  measured  signal  and  noise  surface  wave 

i 

magnitudes,  one  must  have  an  m^  value  for  each  event.  Since  the  parameters 
of  the  events  in  the  data  base  were  originally  derived  from  the  PDE,  NORSAR, 
and  LASA  event  bulletins,  and  since  the  m^  values  of  these  bulletins  often 
disagree,  one  with  the  other,  it  was  necessary  to  find  a means  of  assigning 
an  m^  value  to  each  event  such  that  event-to- event  variance  would  be 
minimized. 

The  first  step  in  assigning  these  m^  values  was  to  form  a 

'combined  m value'  for  each  event.  If  only  one  of  the  event  bulletins  report- 
b 

ed  a given  event,  the  combined  m^  was  simply  the  m^  derived  from  that  bul- 
letin. If  two  event  bulletins  reported  an  event,  the  average  of  the  two  report- 
ed magnitudes  served  as  the  combined  m^.  If  all  three  of  the  lists  reported 
the  event,  the  two  magnitudes  with  the  closest  values  were  averaged  to  form 
the  combined  m^.  Since  the  LASA  event  bulletin  reported  the  greatest  num- 
ber of  events,  the  'combined  m^'  values  are  closer  to  LASA  m^  values  than 
the  PDE  or  NORSAR  m^  values. 

At  this  point,  there  were  four  sets  of  m^  values  assigned  to 

the  events  of  the  data  base:  PDE,  NORSAR,  LASA,  and  combined.  Only  the 

combined  m,  values  had  a magnitude  for  every  event  of  the  data  base.  In 
b 

order  to  decide  which  set  of  values  should  be  used  in  'binning'  the  surface 
wave  magnitudes,  maximum  likelihood  detection  capability  curves  were  com- 
puted for  the  surface  wave  detection  statistics  using  each  of  the  four  types  of 
m^.  The  idea  behind  this  approach  is  that  the  more  precisely  an  m^  estima- 
tion procedure  represents  the  events,  the  less  variance  will  appear  in  the 
maximum  likelihood  fits  to  the  detection  statistics.  The  four  fits  to  the  Ray- 
leigh wave  reference  site  detection  statistics  are  shown  in  Figure  II- 2 to  il- 
lustrate the  approach.  The  upper  portion  of  each  sub-figure  shows  the  detec- 
tion statistics  for  each  data  set.  The  lower  portion  of  each  figure  shows  the 
maximum  likelihood  fit  (solid  line)  to  the  detection  statistics.  The  dashed 
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lines  represent  the  95  percent  confidence  limits  on  this  fit.  This  method  of 
presenting  detection  statistics  is  described  by  Ringdal  (1974).  From  this 
figure  one  can  see  that  the  PDE  and  NORSAR  bulletins  did  not  assign  m^ 
values  to  enough  events  of  the  data  base.  Of  the  four  sets  of  data,  the  com- 
bined m^  values  appear  to  be  the  best  in  terms  of  internal  consistency.  To 
illustrate  this  matter  further,  the  fitting  errors  at  the  50  percent  and  90  per- 
cent detection  thresholds  are  listed  in  Table  II- 1.  The  smallest  errors  occur 


when  the  combined  values  are  used,  indicating  that  the  combined 
best  represent  the  relative  size  of  the  events. 


m values 
b 


Having  measured  the  signal  and  noise  surface  wave  magnitudes 
for  the  reference  site  and  beamformed  array,  and  having  assigned  each  event 
a bodywave  magnitude,  it  is  now  possible  to  test  Ringdal's  method  for  remov- 
ing magnitude  bias. 


SECTION  in 
DATA  ANALYSIS 


A.  REMOVAL  OF  SURFACE  WAVE  MAGNITUDE  BIAS 

The  first  step  in  the  analysis  of  the  data  was  to  run  Ringdal's 
maximum  likelihood  magnitude  estimation  program  using  the  long-period  sur- 
face wave  magnitude  measurements  made  on  detected  signals  and  the  'noise 
magnitude'  measurements  made  on  non-detected  signals.  The  program  oper- 
ates by  binning  these  magnitude  measurements  on  the  basis  of  the  bodywave 
magnitude  of  each  event.  Thus,  for  each  m^  value,  there  is  a set  of  M^ 

values.  Each  set  of  M values  is  processed  by  an  algorithm  designed  to  esti- 

s 

mate  the  maximum  likelihood  magnitude  and  associated  standard  deviation. 

The  algorithm  uses  as  its  starting  point  the  mean  of  the  detected  surface 
waves  and  standard  deviation  set  equal  to  0.  35.  The  algorithm  then  proceeds 
to  determine  the  true  magnitude  and  associated  standard  deviation  by  an  itera- 
tive procedure  which  increments  and  decrements  the  magnitude  and  standard 
deviation  until  the  likelihood  function  (equation  (II-4))  reaches  a maximum 
value.  The  resulting  estimate  of  the  unbiased  surface  wave  magnitude  of  a 
set  of  events  observed  at  one  station  is  operationally  equivalent  to  estimating 
the  surface  wave  magnitude  of  a single  event  from  station  magnitude  measure- 
ments made  on  a network  of  seismic  stations. 

The  output  of  this  program  is  plotted  in  Figure  HI- 1 for  vertical 
Rayleigh  and  in  Figure  QI-Z  for  transverse  Love  surface  wave  measurements. 
The  open  circles  in  these  figures  represent  the  reference  site  maximum  like- 
lihood estimate  of  the  surface  wave  magnitude.  The  solid  dots  represent  the 
corresponding  reference  site  mean  signal  surface  wave  magnitudes.  The 
crosses  represent  the  beamformed  array  mean  signal  surface  wave  magnitudes 
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which  are  to  serve  as  the  standard  against  which  to  judge  whether  the  refer- 
ence site  mean  signal  magnitudes  contain  positive  magnitude  bias  and,  if  so, 
whether  Ringdal's  method  removes  this  bias.  Finally,  the  triangles  repre- 
sent the  beamformed  array  maximum  likelihood  estimate  of.  the  surface  wave 
magnitude. 

Examination  of  these  figures  indicates: 

• As  a result  of  averaging  only  the  surface  wa/e  magnitudes  of 
detected  events  (forming  conventional  network  estimates), 
significant  positive  bias  is  observed. 

• This  bias  is  significant  below  the  reference  site  (single-sensor) 
90  percent  detection  threshold. 

• Based  on  a comparison  of  ALPA  single- sensor  estimates 

to  ALPA  array  M estimates,  Ringdal's  unbiased  magnitude 

8 

estimation  method  accurately  corrects  single- sensor  magni- 
tudes down  to  approximately  the  25  percent  detection  threshold. 

• There  are  indications  of  overcorrection  by  Ringdal's  method 

below  the  25  percent  detection  threshold,  resulting  in  under- 
estimation of  The  possibility  of  this  effect  was  indicated 

by  von  Seggern  and  Rivera  (1978). 

• The  observation  of  valid  bias  removal  by  Ringdal's  method  in- 
dicates that  his  assumption  that  magnitude  deviations  occur 
independently  of  noise  level  deviations  is  applicable  to  estimat- 
ing unbiased  mean  surface  wave  magnitudes  corresponding  to 
m^  values. 

The  reference  site  mean  signal  and  maximum  likelihood  esti- 
mate values  at  m^=  4.  2 and  at  m^  values  below  the  reference  site  25  per- 
cent detection  threshold  show  large  departures  from  the  trend  of  the  surround- 
ing data.  For  these  points,  there  were  only  one  or  two  detections  and  ten  to 
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twenty  non- detections  in  the  event  population.  This  indicates  that  Ringdal's 
method  may  not  work  well  when  the  detection  population  is  sparse.  There- 
fore, these  data  points  will  be  dropped  from  the  remainder  of  the  analysis. 

For  m,  values  where  there  were  no  surface  wave  non-detections, 
b 

the  beamformed  array  mean  signal  M values  trend  on  the  average  approxi- 

8 

mately  0.  08  M units  below  the  corresponding  reference  site  M values.  This 
s s 

is  to  be  expected,  since  the  beamforming  process  attenuates  the  signal  ampli- 
tudes slightly.  In  the  evaluation  of  ALPA  (Strauss,  1973)  this  signal  attenua- 
tion was  estimated  to  be  approximately  1 dB  or  0.  05  magnitude  units. 

From  Figures  HI- 1 and  III- 2,  it  appears  that  Ringdal's  method 
does  remove  positive  magnitude  bias  from  surface  wave  magnitude  data.  How- 
ever, each  estimate  shown  in  the  two  figures  carries  with  it  a standard 
deviation.  (Examples  of  these  are  listed  in  Table  IH-1  for  the  vertical  Ray- 
leigh data.  ) If  these  standard  deviations  are  very  large,  one  could  not  state 
with  any  degree  of  certainty  that  the  mean  signal  magnitudes  were  not  equal 
to  the  maximum  likelihood  estimate  magnitudes.  To  resolve  this  matter,  a 
hypothesis  test  procedure  was  carried  out. 

Following  the  procedure  detailed  by  Johnson  and  Leone  (1964), 
the  hypothesis  that  the  reference  site  mean  signal  Mg  equals  the  reference 
site  maximum  likelihood  was  tested  using  the  Student's  t statistic.  The 
test  statistic  t was  computed  from 

(ni-1) 


with  v degrees  of  freedom,  where 
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TABLE  m-1 

SURFACE  WAVE  MAGNITUDES  WITH  ASSOCIATED 
STANDARD  DEVIATIONS 
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nl  + 1 


n2  + i 


In  these  equations 


x.  = M estimate 
1 s 

S.  = sample  standard  deviation  of  x. 
i l 

n^  = number  of  points  used  to  compute  and  x.  . 

The  results  of  this  test  showed  that  the  hypothesis  failed,  i.  e.  , 

the  reference  site  mean  signal  M and  the  reference  site  maximum  likelihood 

s 

estimate  were  not  equal,  with  90  percent  confidence,  for  m^  4.  1, 

4.  3,  and  4.  4.  (Recall  that  the  m = 4.  2 data  were  dropped  from  the  analysis.  ) 

Based  on  this  test,  at  m^  values  greater  than  4.4,  the  reference  site  mean 

signal  M and  unbiased  M could  not  be  said  to  be  unequal. 

3 8 

Since  the  beamformed  array  mean  signal  M values  are  the 

8 

standard  of  comparison,  it  is  important  to  perform  this  test  with  these  values 

and  the  reference  site  maximum  likelihood  estimate  M values.  In  this  case, 

s 

the  computed  values  of  the  t statistic  were  such  that  at  no  m^  value  could 
these  two  estimates  of  be  said  to  be  unequal. 

The  above  analysis  indicates  that  the  conventional  Mg  estimates 
(average  of  the  individual  signal  values)  do  contain  a positive  bias  for  Mg 
measured  on  events  with  m^  below  the  90  percent  detection  threshold  and  that 
Ringdal's  method  does  remove  this  bias,  at  least  down  to  the  beamformed  ar- 
ray 50  percent  detection  threshold  (equivalent  to  the  reference  site  25  percent 
detection  threshold).  This  implies  an  extension  of  the  Mg-  m^  curve  by  about 
0.  8 m^  units  for  this  data. 


HI-7 


The  statistical  tests  indicate  that  the  observed  bias  corrections 
were  significant  at  a 90  percent  confidence  level.  No  significant  difference 
could  be  observed  between  the  bias  corrected  single  sensor  reference  site 
magnitudes  and  the  presumed  unbiased  array  magnitude  measurements,  also 
at  a confidence  level  of  90  percent. 

It  is  of  interest  to  examine  the  least- squares  line  fits  to  the 

various  data  sets.  In  each  of  the  following  figures,  one  line  was  fit  through 

the  m^  points  from  = 4.  0 to  m^  = 4.  9 and  a second  line  through  the 

m^  points  from  m^  = 4.  9 to  = 6.  2.  These  fits  were  required  to  tie 

at  the  intersecting  point.  The  selection  of  the  m^=  4.  9 intersection  point  was 

guided  by  the  observation  that  positive  magnitude  bias  was  significant  below 

m^=  4.  9 for  reference  site  data.  Thus,  breaking  the  Mg-  mfa  data  into  two 

sets  at  this  point  should  most  clearly  demonstrate  the  effect  of  this  bias  on 

the  M^-  m^  relationship.  Fits  to  the  beamformed  array  M - m^  data  were 

made  over  the  same  ranges  in  m^  for  ease  of  comparison  with  the  reference 

site  M - m,  fits, 
s b 

In  the  fitting  procedure,  each  mean  or  unbiased  M - m^  point 

was  weighted  by  the  number  of  individual  Mg  values  used  to  compute  the  mean 

or  unbiased  M value,  The  standard  deviation  associated  with  each  mean  or 
8 

unbiased  value  was  likewise  weighted  to  produce  the  standard  deviation  of 
of  the  Mg-  m^  fit.  The  standard  deviations  so  computed  are  larger  than 
one  would  expect  from  the  Mg-  m^  points  shown  in  the  following  figures, 
since  the  Mg-  m^  points  shown  are  the  mean  values  representing  the  M pop- 
ulation at  each  m^  value. 


Figure  III- 3a  shows  the  reference  site  mean  signal  M - 


mb 


data,  while  Figure  HI- 3b  shows  the  corresponding  unbiased  M - m^  data. 
Figure  III-4a  shows  the  beamformed  array  mean  signal  M#-  m^  data,  while 
Figure  ni-4b  shows  the  corresponding  unbiased  M(-  m^  data.  Although  not 
used  in  the  fitting  procedure,  the  points  below  m^=  4.  0 are  shown  for  the 
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sake  of  completeness.  These  figures  all  show  the  vertical  component  20- 
second  Rayleigh  wave  Mfl-  m^  data.  The  corresponding  figures  for  transverse 
component  Love  wave  M^-  m^  data  are  not  shown  since  they  are  essentially 
identical  to  the  Rayleigh  wave  data.  The  equations  for  the  Mg-  m^  fits  shown 
in  Figures  III- 3 and  HI-4  are  listed  in  Table  III-2  along  with  their  associated 
standard  deviations.  Note  that  these  standard  deviations  are  all  essentially 
the  same,  with  the  exception  of  the  standard  deviation  (0.  27)  associated  with 
equation  (HI- 3).  Since  this  equation  is  based  on  the  smallest  amount  of  data, 
it  is  felt  that  this  difference  is  not  significant.  For  comparison,  note  that 
Lambert  et  al.  (1974)  found  a value  of  0.  37  for  the  standard  deviation  of  fits 
to  HGLP  m^  data,  while  Strauss  and  Weltman  (1977)  found  a value  of 
0.  35  for  the  standard  deviation  of  fits  to  SRO  Mg-  data. 

Referring  to  Figure  HI- 3a,  note  that  the  data  show  a pronounced 
change  of  slope  at  m^  4.  9,  indicating  magnitude  bias  below  this  point  for  the 
reference  site  mean  signal  Mg-  m^  data.  The  corresponding  unbiased  M -m^ 
data  of  Figure  HI- 3b  shows  a much  smaller  change  in  slope.  The  beamform- 
ed  array  mean  signal  Mg-  m^  data  of  Figure  HI-4a  show  a small  change  in 
slope,  indicating  the  potential  presence  of  a small  amount  of  magnitude  bias. 
Again,  the  corresponding  maximum  likelihood  beam  estimate  Mg-  m^  data  of 
Figure  HI-4b  show  a smaller  change  in  slope.  In  fact,  comparison  of  Figures 
HI- 3b  and  HI-4b  (which  shows  essentially  parallel  Mg-  m^  curves),  indicates 
that  Ringdal's  method  corrects  the  reference  site  Mg  data  for  positive  magni- 
tude bias  for  mfe  values  between  the  90  percent  and  the  25  percent  detection 
thresholds.  This  results  in  an  extension  of  the  linear  Mg-  m^  relationship 
from  approximately  the  reference  site  90  percent  detection  threshold  to  the 
reference  site  25  percent  detection  threshold.  For  this  case,  it  is  shown 
that  an  unbiased  reference  curve  could  be  represented  down  to  4.  0 using 
Ringdal's  method  versus  4.  8 by  averaging  the  of  detected  long-period 
signals. 
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One  further  point  remains  to  be  made  concerning  these  Mg-m^ 

fits.  In  the  past  it  has  been  common  practice  to  compute  a linear  fit  to  the 

Mg-  m^  data  throughout  the  range  in  m^  of  the  data  used.  This  results  in 

m^  relationships  which  are  dependent  on  this  range  in  m^.  For  example, 

consider  the  following  M - m relationships: 

s b 


M = 1.18 
8 

M =1.36 


- 1.  66 


1. 44 


(Lambert  et  al.  , 1974)  (III- 11) 


1.  59  nr  - 3.  97 
b 


(Basham,  1969) 
(Richter,  1958)  . 
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Equation  (HI-11)  was  derived  from  ALPA,  NORSAR,  and  HGLP 
surface  wave  data  with  a center  of  mass  at  m = 4.  37.  This  equation  is  almost 
identical  to  equation  (HI- 9)  of  Table  IH-2,  which  is  defined  over  the  range 
4.  0 ^ m^  < 4.  9.  Equation  (HI-12),  with  a center  of  mass  at  m^3  4.  85  and  a 
slope  midway  between  those  of  equations  (HI- 9)  and  (HI- 10)  of  Table  IH-2, 
may  be  considered  to  be  essentially  an  average  fit  to  the  data  used  to  compute 
equations  (HI-9)  and  (HI-10).  Equation  (HI-13),  which  was  computed  from 
large  m events  (m  > 6.  5),  is  almost  identical  to  equation  (HI-10),  which  was 
computed  for  the  range  4.  9 ^ <6.2.  All  of  this  indicates  that  the  M^-  m^ 

data  is  best  fitted  by  two  lines.  While  the  evidence  is  not  conclusive,  it  ap- 
pears that  the  point  of  intersection  of  these  fits  should  be  at  about  m^s  5.  0. 

Although  the  results  of  the  unbiased  estimates  of  Mg- 
relationship  are  consistent  with  other  authors'  results,  some  idea  of  the 
physical  significance  and  generality  of  the  results  can  be  obtained  by  com- 
parison with  bodywave  source  spectrum  studies.  In  a study  of  P wave  spec- 
tral measurements  by  Sax  (1975),  the  low  frequency  spectral  displacement 
amplitude,  defined  as  the  spectrum  below  the  corner  frequency,  scales  with 
m.  as  follows. 


Log  S3  = 1.14  m,  - 1.  21 


(HI- 14) 


HI- 13 


where  Log  is  the  spectral  magnitude  corrected  for  distance  as 


- • 


Log  = Log  Qq  + k log  A + C (HI-15) 


where 


regional  k = 3,  C = -1.  59 
teleseismic  k = 1,  C = 1.89. 

Log  is  the  measured  log  (displacement  amplitude),  and  A is  the  epicentral 
distance  of  propagation.  Aki  (1967)  treats  M measurement  as  equivalent  to 
the  spectral  magnitude  of  surface  waves  to  within  an  additive  constant.  Com- 
paring equation  (HI-14)  with  equation  (IH-9),  it  is  observed  that  the  low  fre- 
quency spectral  magnitude  of  bodywaves  also  scales  similarly  to  M to  within 
a constant,  0.  58,  or  approximately  a factor  of  4.  This  can  be  shown  to  be 
reasonable  by  observing  the  scaling  relationship  for  computing  magnitude 
from  spectrum,  given  by  Thatcher  and  Hanks  (1973)  for  local  magnitude  by 
defining  local  magnitude  as  the  square  root  of  the  propagated  energy. 


M = Log  Q + - log  f . 
L o 2 c 


(in- 16) 


Sax  (1975)  extended  the  spectral  scaling  law  to  the  magnitude  of  P waves, 
m^  , and  teleseismic  distances  as 


mbs  " ***  °o  + 2 108  fc 


(HI- 17) 


by  means  of  the  attenuation  function  correction  given  by  equation  (HI- 15), 
which  was  determined  by  fitting  data.  Sax  (1975)  observed  that  the  average 
logarithmic  corner  frequency  of  P waves  decreases  with  m^  as 


Logfc  = 


•0. 074 


“b 


+ 0.  74 


(HI- 18) 
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which  leads  to 


Log  Qo  = 1.11  -1.11  . (HI-19) 

The  predicted  scaling  of  the  low  frequency  spectral  magnitude  of  P waves  is 
reasonably  close  to  the  observed  scaling  law  given  by  equation  (HI- 14).  Thus, 
the  small  increment  in  the  slope  of  the  M^-  m^  relationship  probably  reflect* 
the  decrease  in  the  comer  frequencies  of  larger  magnitude  events. 

The  results  for  scaling  bodywave  spectral  amplitudes  did  differ 
from  the  results  in  the  Mfl-  relationship  in  that  no  break  in  the  slope  was 
observed.  This  might  be  explained  by  insufficient  observations  of  events 
above  rn^=  5 to  determine  the  break  in  slope  of  bodywave  spectral  amplitude 
versus  m^,  should  it  exist. 

The  break  in  slope  of  the  Mg-  m^  relationship  may  be  due  to 
the  effect  of  the  short-period  system  response  curve  on  the  measurement  of 
m^.  The  corrections  for  system  magnification  would  be  expected  to  be  more 
accurate  for  smaller  events  with  corner  frequencies  near  the  peak  of  the  sys- 
tem response.  Errors  in  magnification  corrections  could  be  introduced  pro- 
gressively in  larger  magnitude  events  due  to  bias  in  visual  measurement  of 
period  of  a signal.  The  bias  would  be  toward  measuring  shorter  periods  than 
would  be  the  case  if  the  ground  displacement  wavelets  were  observed  directly. 
This  could  result  in  compression  of  the  scale  above  00^=  5. 

B.  THE  EFFECT  OF  MAGNITUDE  BIAS  CORRECTION  ON  THE  M - rm 
DISCRIMINANT  8 

Figure  IH-5  shows  a typical  Mg-  m^  scatter  plot.  The  stars 
represent  the  M#-  m^  values  for  the  Kurile  Islands  - Kamchatka  events  of  the 
data  base  as  measured  on  the  beamformed  array  traces.  The  open  circles 
represent  Mg-  values  for  a suite  of  Eurasian  presumed  nuclear  explosions 
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their  m^  trend  in  parallel  with  the  earthquake  data  for  at  least  one-half 
magnitude  unit.  It  is  believed  that  the  non-parallelism  of  the  lines  fit  to  the 
upper  m^  range  of  earthquake  and  presumed  explosion  data  is  probably  esti- 
mation error  due  to  the  very  small  size  of  the  available  presumed  explosion 
population. 

Since  it  will  always  be  difficult  to  establish  a large  presumed 
nuclear  explosion  data  base,  an  alternative  to  the  preceding  plots  is  offered 
as  a means  of  utilizing  this  discriminant.  This  approach  consists  of  first 
establishing  network  Mfl  values  for  a large  number  of  earthquakes  to  establish 
the  earthquake  Mg-  m^  relationship  for  the  network.  A test  event  can  be  com- 
pared to  this  earthquake  population  by  computing  the  separation  S between 
the  earthquake  network  derived  from  the  established  M - m^  relationship 
and  the  measured  network  M of  the  test  event.  If  this  separation  is  greater 
than  some  threshold  (for  example,  greater  than  two  standard  deviations  of  the 
earthquake  M^-  m^  fit),  the  test  event  can  be  identified  as  a presumed  nuclear 
explosion. 

This  procedure  was  carried  out  for  the  mean  signal  data 

and  the  maximum  likelihood  estimate  M data  computed  from  the  vertical 

8 

Rayleigh  wave  measurements  made  on  beamformed  array  data.  The  M#-  m^ 
relationships  used  were  equations  (HI- 7)  and  (III- 8)  for  the  mean  signal 
data  and  equations  (HI-9)  and  (HI-10)  for  the  maximum  likelihood  estimate 
data.  The  value  of  two  standard  deviations  was  set  as  the  discrimination 
threshold  because  95  percent  of  the  members  of  a normally  distributed  popula- 
tion lie  within  two  standard  deviations  of  the  mean.  Thus,  the  probability  that 
a test  event  is  a member  of  the  earthquake  population  when  it  is  separated  from 
the  earthquake  mean  by  a positive  deviation  of  more  than  two  standard  devia- 
tions is  less  than  0.  025.  In  this  example,  a value  of  0.  4 was  used  as  the 
standard  deviation. 
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In  order  to  provide  more  ineight  into  the  effect  on  the  Mj- 
discriminant  of  removing  magnitude  bias,  the  separation  S is  presented  in 
Table  III- 3 for  the  mean  signal  Mg  (S  . ) and  the  maximum  likelihood  estimate 
(Smje).  This  table  shows  that  for  a discrimination  threshold  of  two  stand- 
ard deviations: 


• In  neither  case  are  any  of  the  earthquakes  mis -identified. 

• Eleven  of  the  fifteen  presumed  nuclear  explosions  are  so 
identified  when  the  mean  signal  M is  used. 

• Twelve  of  the  fifteen  presumed  nuclear  explosions  are  so 

identified  when  the  unbiased  M is  used. 

s 

• Those  presumed  explosions  which  were  not  identified  as  such 
are  the  smaller  m^  events  for  which  there  was  only  one  obser- 
vation available.  The  one  event  which  was  identified  as  an 
earthquake  using  mean  signal  M and  as  an  explosion  using  un- 
biased  was  one  of  the  three  for  which  there  were  sufficient 
observations  to  apply  Ringdal's  method.  (The  other  two  of  the 

three  were  identified  as  explosions  using  both  types  of  M . ) 

s 

If  a more  conservative  threshold  of  three  standard  deviations 
is  used,  only  seven  of  the  fifteen  presumed  nuclear  explosions  are  identified 
when  the  mean  signal  is  used,  while  ten  are  identified  as  explosions  when 

the  unbiased  M is  used. 

0 

These  results  indicate  that  Ringdal's  method  of  removing  posi- 
tive magnitude  bias  from  network  m^  can  improve  the  value  of  the  M - 
discriminant.  Determination  of  the  degree  by  which  it  will  improve  the  dis- 
criminant must  await  the  collection  of  a presumed  explosion  data  base  con- 
taining multiple  observations  of  each  event,  i.  e. , by  a network  of  seismic 
stations. 
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TABLE  in- 3 

SEPARATION  S FROM  EARTHQUAKE  MEAN  M 


Earthquakes 

Presumed  Nuclear  Explosions 

""b 

S . 

Big 

Smle 

"'b 

S . 

Smle 

0.  02 
0.00 
0.  05 
-0.  11 
0.09 
-0.12 
-0.  26 
0.  05 
-0.  07 


0. 14 
0.08 
-0.  02 
0.16 
0.  06 
-0.43 
0.  45 
-0.  22 
-0.  02 


0.05 
0.  02 
-0.  01 
-0. 13 
0.  00 
-0.  20 
-0.  16 
0.02 
-0.08 
0.  14 
0.  14 
0.  08 
-0.  02 
0.  16 
0.  06 
-0.44 
0.44 
-0.23 
-0.03 


0.  39 
0.  83 


0.  87 
1. 30 

0. 69 
1. 11 
1.43 
1.46 

1.  57 
1.49 
1.51 
1.  89 


1.  30 

0.  83 
1.13 
1. 74 

1.  29 

0.  69 

1.  82 
1.42 
1.45 
1.  56 
1.48 
1.  50 
1.  88 


S . denotes  separation  between  mean  signal  Ms-m^  values  and  the 
® line  fit  to  the  earthquake  mean  signal  Mg-m^  data 

Smie  denotes  separation  between  unbiased  Ms-m^  values  and  the  line 
fit  to  the  earthquake  unbiased  Mg-m^  data 

One  standard  deviation  is  taken  as  0.  4 magnitude  units. 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  following  conclusions  have  been  reached  in  the  course  of 
this  test  of  Ringdal's  method  of  removing  positive  magnitude  bias: 

• Surface  wave  magnitudes  are  strongly  affected  by  positive 
magnitude  bias  for  events  with  bodywave  magnitudes  below  the 
90  percent  detection  threshold. 

• Below  this  detection  threshold,  this  bias  is  the  controlling  fac- 
tor in  the  trend  of  the  Mg-  m^  relationship. 

• Ringdal's  method  of  removing  this  bias  successfully  removed 
the  bias  from  the  surface  wave  magnitudes  for  events  with 
bodywave  magnitudes  down  to  approximately  the  25  percent 
detection  threshold. 

• For  good  results,  the  method  requires  at  least  three  detections 
with  corresponding  signal  magnitudes. 

• The  earthquake  M#-  m^  relationships  derived  from  unbiased 

ALPA- recorded  surface  wave  magnitude  data  are: 

(M  s 1.09  m.  - 1.47  4.  0 < m<  4.  9 

Single  Site  / 8 b 

(M  =1.53  m_  - 3.  61  4.9£m.<6.2 


Array 


!M  =1.14 

* 

M = 1. 


- 1.79  4.0Sm,  <4.9 


- 3.  79  4.  9<  m,  < 6.  2 . 


• Removal  of  magnitude  bias  from  those  presumed  explosions 

for  which  there  were  multiple  observations  brought  the  M(  of 

these  events  into  line  with  the  trend  of  the  M data  for  those 

s 

presumed  explosions  with  m^  at  or  above  the  array  100  per- 
cent detection  threshold  of  surface  waves. 


B. 


RECOMMENDATIONS 


In  light  of  the  preceding  discussion  of  positive  magnitude  bias 
removal  and  Mg-  m^  relationships,  and  considering  the  difficulty  of  assembling 
a presumed  nuclear  explosion  data  base  large  enough  to  establish  a reliable 
explosion  M^-  m^  relationship,  it  would  appear  that  the  optimum  manner  in 
which  to  use  the  Mg-  m^  discriminant  would  be  the  determination  of  the  sep- 
aration between  each  of  the  M^-  points  and  the  earthquake  M(-  m^  relation- 
ship as  described  in  Section  III.  This  approach  permits  one  to  assign  a single 
number  as  the  discrimination  value  rather  than  a pair  of  numbers  (Mg  and  m^). 
This  discrimination  number  can  then  be  compared  to  some  pre-set  threshold 
(such  as  two  or  three  standard  deviations  from  the  mean)  in  order  to  identify 
the  event. 


It  is  recommended  that  a test  of  this  procedure  using  unbiased 
Mg  data  be  carried  out  using  the  SRO-ASRO  stations  as  a network.  From  the 
signal  and  noise  magnitudes  measured  at  the  stations  for  a suite  of  earthquakes, 
the  unbiased  Mg-  m^  relationship  can  be  established.  Then  a suite  of  test 
events  can  be  used  to  determine  the  discrimination  capability  of  this  method 
in  terms  of  the  percentage  of  correct  event  identifications  and  the  m^  range 
over  which  it  is  valid. 

The  results  of  this  study  also  indicate  that  an  in-depth  investiga- 
tion of  M#-  m^  relationships  using  unbiased  mean  M#  data  would  be  fruitful. 

In  particular,  more  large  (m^  > 5.  0)  earthquakes  are  needed  to  better  estab- 
lish the  M(-  mfe  least- squares  trend  line  for  large  m^.  The  key  points  to  be 
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established  are  whether  Mg-  data  is  best  fitted  with  two  lines,  as  illustrated 
| in  this  report  and,  if  so,  where  the  point  of  intersection  lies.  The  same  points 

need  to  be  answered  for  explosion  Mg-  m^  data.  The  difficulty  of  assembling 
| a sufficiently  large  explosion  data  base  at  any  one  station  could  be  circum- 

vented by  removing  the  station  instrument  response  and  by  comparing  similar 
| measurements  of  Mg  at  all  stations  at  all  times.  A study  of  this  nature  can 

be  used  to  test  the  hypothesis  that  the  earthquake  and  explosion  M - m^  rela- 
| tionships  are  parallel. 

I 
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